Proteome-wide base editor screens to assess phosphorylation site functionality in high-throughput.

Signaling pathways that drive gene expression are typically depicted as having a dozen or so landmark phosphorylation and transcriptional events. In reality, thousands of dynamic post-translational modifications (PTMs) orchestrate nearly every cellular function, and we lack technologies to find causal links between these vast biochemical pathways and genetic circuits at scale. Here, we describe "signaling-to-transcription network" mapping through the development of PTM-centric base editing coupled to phenotypic screens, directed by temporally-resolved phosphoproteomics. Using T cell activation as a model, we observe hundreds of unstudied phosphorylation sites that modulate NFAT transcriptional activity. We identify the phosphorylation-mediated nuclear localization of the phosphatase PHLPP1 which promotes NFAT but inhibits NFκB activity. We also find that specific phosphosite mutants can alter gene expression in subtle yet distinct patterns, demonstrating the potential for fine-tuning transcriptional responses. Overall, base editor screening of PTM sites provides a powerful platform to dissect PTM function within signaling pathways.

Optimization of Cas12a for multiplexed genome-scale transcriptional activation.

Cas12a CRISPR technology, unlike Cas9, allows for facile multiplexing of guide RNAs from a single transcript, simplifying combinatorial perturbations. While Cas12a has been implemented for multiplexed knockout genetic screens, it has yet to be optimized for CRISPR activation (CRISPRa) screens in human cells. Here, we develop a new Cas12a-based transactivation domain (TAD) recruitment system using the ALFA nanobody and demonstrate simultaneous activation of up to four genes. We screen a genome-wide library to identify modulators of growth and MEK inhibition, and we compare these results with those obtained with open reading frame (ORF) overexpression and Cas9-based CRISPRa. We find that the activity of multiplexed arrays is largely predictable from the best-performing guide and provide criteria for selecting active guides. We anticipate that these results will greatly accelerate the exploration of gene function and combinatorial phenotypes at scale.

Scanning mutagenesis of the voltage-gated sodium channel NaV1.2 using base editing.

It is challenging to apply traditional mutational scanning to voltage-gated sodium channels (NaVs) and functionally annotate the large number of coding variants in these genes. Using a cytosine base editor and a pooled viability assay, we screen a library of 368 guide RNAs (gRNAs) tiling NaV1.2 to identify more than 100 gRNAs that change NaV1.2 function. We sequence base edits made by a subset of these gRNAs to confirm specific variants that drive changes in channel function. Electrophysiological characterization of these channel variants validates the screen results and provides functional mechanisms of channel perturbation. Most of the changes caused by these gRNAs are classifiable as loss of function along with two missense mutations that lead to gain of function in NaV1.2 channels. This two-tiered strategy to functionally characterize ion channel protein variants at scale identifies a large set of loss-of-function mutations in NaV1.2.

Chromatin complex dependencies reveal targeting opportunities in leukemia.

Chromatin regulators are frequently mutated in human cancer and are attractive drug targets. They include diverse proteins that share functional domains and assemble into related multi-subunit complexes. To investigate functional relationships among these regulators, here we apply combinatorial CRISPR knockouts (KOs) to test over 35,000 gene-gene pairings in leukemia cells, using a library of over 300,000 constructs. Top pairs that demonstrate either compensatory non-lethal interactions or synergistic lethality enrich for paralogs and targets that occupy the same protein complex. The screen highlights protein complex dependencies not apparent in single KO screens, for example MCM histone exchange, the nucleosome remodeling and deacetylase (NuRD) complex, and HBO1 (KAT7) complex. We explore two approaches to NuRD complex inactivation. Paralog and non-paralog combinations of the KAT7 complex emerge as synergistic lethal and specifically nominate the ING5 PHD domain as a potential therapeutic target when paired with other KAT7 complex member losses. These findings highlight the power of combinatorial screening to provide mechanistic insight and identify therapeutic targets within redundant networks.

Accounting for small variations in the tracrRNA sequence improves sgRNA activity predictions for CRISPR screening.

CRISPR technology is a powerful tool for studying genome function. To aid in picking sgRNAs that have maximal efficacy against a target of interest from many possible options, several groups have developed models that predict sgRNA on-target activity. Although multiple tracrRNA variants are commonly used for screening, no existing models account for this feature when nominating sgRNAs. Here we develop an on-target model, Rule Set 3, that makes optimal predictions for multiple tracrRNA variants. We validate Rule Set 3 on a new dataset of sgRNAs tiling essential and non-essential genes, demonstrating substantial improvement over prior prediction models. By analyzing the differences in sgRNA activity between tracrRNA variants, we show that Pol III transcription termination is a strong determinant of sgRNA activity. We expect these results to improve the performance of CRISPR screening and inform future research on tracrRNA engineering and sgRNA modeling.

Bidirectional genome-wide CRISPR screens reveal host factors regulating SARS-CoV-2, MERS-CoV and seasonal HCoVs.

CRISPR knockout (KO) screens have identified host factors regulating severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) replication. Here, we conducted a meta-analysis of these screens, which showed a high level of cell-type specificity of the identified hits, highlighting the necessity of additional models to uncover the full landscape of host factors. Thus, we performed genome-wide KO and activation screens in Calu-3 lung cells and KO screens in Caco-2 colorectal cells, followed by secondary screens in four human cell lines. This revealed host-dependency factors, including AP1G1 adaptin and ATP8B1 flippase, as well as inhibitors, including mucins. Interestingly, some of the identified genes also modulate Middle East respiratory syndrome coronavirus (MERS-CoV) and seasonal human coronavirus (HCoV) (HCoV-NL63 and HCoV-229E) replication. Moreover, most genes had an impact on viral entry, with AP1G1 likely regulating TMPRSS2 activity at the plasma membrane. These results demonstrate the value of multiple cell models and perturbational modalities for understanding SARS-CoV-2 replication and provide a list of potential targets for therapeutic interventions.

Benchmarking of SpCas9 variants enables deeper base editor screens of BRCA1 and BCL2.

Numerous rationally-designed and directed-evolution variants of SpCas9 have been reported to expand the utility of CRISPR technology. Here, we assess the activity and specificity of WT-Cas9 and 10 SpCas9 variants by benchmarking their PAM preferences, on-target activity, and off-target susceptibility in cell culture assays with thousands of guides targeting endogenous genes. To enhance the coverage and thus utility of base editing screens, we demonstrate that the SpCas9-NG and SpG variants are compatible with both A > G and C > T base editors, more than tripling the number of guides and assayable residues. We demonstrate the performance of these technologies by screening for loss-of-function mutations in BRCA1 and Venetoclax-resistant mutations in BCL2, identifying both known and new mutations that alter function. We anticipate that the tools and methodologies described here will facilitate the investigation of genetic variants at a finer and deeper resolution for any locus of interest.

Bidirectional genome-wide CRISPR screens reveal host factors regulating SARS-CoV-2, MERS-CoV and seasonal HCoVs.

Several genome-wide CRISPR knockout screens have been conducted to identify host factors regulating SARS-CoV-2 replication, but the models used have often relied on overexpression of ACE2 receptor. Additionally, such screens have yet to identify the protease TMPRSS2, known to be important for viral entry at the plasma membrane. Here, we conducted a meta-analysis of these screens and showed a high level of cell-type specificity of the identified hits, arguing for the necessity of additional models to uncover the full landscape of SARS-CoV-2 host factors. We performed genome-wide knockout and activation CRISPR screens in Calu-3 lung epithelial cells, as well as knockout screens in Caco-2 intestinal cells. In addition to identifying ACE2 and TMPRSS2 as top hits, our study reveals a series of so far unidentified and critical host-dependency factors, including the Adaptins AP1G1 and AP1B1 and the flippase ATP8B1. Moreover, new anti-SARS-CoV-2 proteins with potent activity, including several membrane-associated Mucins, IL6R, and CD44 were identified. We further observed that these genes mostly acted at the critical step of viral entry, with the notable exception of ATP8B1, the knockout of which prevented late stages of viral replication. Exploring the pro- and anti-viral breadth of these genes using highly pathogenic MERS-CoV, seasonal HCoV-NL63 and -229E and influenza A orthomyxovirus, we reveal that some genes such as AP1G1 and ATP8B1 are general coronavirus cofactors. In contrast, Mucins recapitulated their known role as a general antiviral defense mechanism. These results demonstrate the value of considering multiple cell models and perturbational modalities for understanding SARS-CoV-2 replication and provide a list of potential new targets for therapeutic interventions.

Bidirectional genome-wide CRISPR screens reveal host factors regulating SARS-CoV-2, MERS-CoV and seasonal coronaviruses.

Several genome-wide CRISPR knockout screens have been conducted to identify host factors regulating SARS-CoV-2 replication, but the models used have often relied on overexpression of ACE2 receptor. Additionally, such screens have yet to identify the protease TMPRSS2, known to be important for viral entry at the plasma membrane. Here, we conducted a meta-analysis of these screens and showed a high level of cell-type specificity of the identified hits, arguing for the necessity of additional models to uncover the full landscape of SARS-CoV-2 host factors. We performed genome-wide knockout and activation CRISPR screens in Calu-3 lung epithelial cells, as well as knockout screens in Caco-2 intestinal cells. In addition to identifying ACE2 and TMPRSS2 as top hits, our study reveals a series of so far unidentified and critical host-dependency factors, including the Adaptins AP1G1 and AP1B1 and the flippase ATP8B1. Moreover, new anti-SARS-CoV-2 proteins with potent activity, including several membrane-associated Mucins, IL6R, and CD44 were identified. We further observed that these genes mostly acted at the critical step of viral entry, with the notable exception of ATP8B1, the knockout of which prevented late stages of viral replication. Exploring the pro- and anti-viral breadth of these genes using highly pathogenic MERS-CoV, seasonal HCoV-NL63 and -229E and influenza A orthomyxovirus, we reveal that some genes such as AP1G1 and ATP8B1 are general coronavirus cofactors. In contrast, Mucins recapitulated their known role as a general antiviral defense mechanism. These results demonstrate the value of considering multiple cell models and perturbational modalities for understanding SARS-CoV-2 replication and provide a list of potential new targets for therapeutic interventions.

Massively parallel assessment of human variants with base editor screens.

Understanding the functional consequences of single-nucleotide variants is critical to uncovering the genetic underpinnings of diseases, but technologies to characterize variants are limiting. Here, we leverage CRISPR-Cas9 cytosine base editors in pooled screens to scalably assay variants at endogenous loci in mammalian cells. We benchmark the performance of base editors in positive and negative selection screens, identifying known loss-of-function mutations in BRCA1 and BRCA2 with high precision. To demonstrate the utility of base editor screens to probe small molecule-protein interactions, we screen against BH3 mimetics and PARP inhibitors, identifying point mutations that confer drug sensitivity or resistance. We also create a library of single guide RNAs (sgRNAs) predicted to generate 52,034 ClinVar variants in 3,584 genes and conduct screens in the presence of cellular stressors, identifying loss-of-function variants in numerous DNA damage repair genes. We anticipate that this screening approach will be broadly useful to readily and scalably functionalize genetic variants.

Optimization of AsCas12a for combinatorial genetic screens in human cells.

Cas12a RNA-guided endonucleases are promising tools for multiplexed genetic perturbations because they can process multiple guide RNAs expressed as a single transcript, and subsequently cleave target DNA. However, their widespread adoption has lagged behind Cas9-based strategies due to low activity and the lack of a well-validated pooled screening toolkit. In the present study, we describe the optimization of enhanced Cas12a from Acidaminococcus (enAsCas12a) for pooled, combinatorial genetic screens in human cells. By assaying the activity of thousands of guides, we refine on-target design rules and develop a comprehensive set of off-target rules to predict and exclude promiscuous guides. We also identify 38 direct repeat variants that can substitute for the wild-type sequence. We validate our optimized AsCas12a toolkit by screening for synthetic lethalities in OVCAR8 and A375 cancer cells, discovering an interaction between MARCH5 and WSB2. Finally, we show that enAsCas12a delivers similar performance to Cas9 in genome-wide dropout screens but at greatly reduced library size, which will facilitate screens in challenging models.

Genome-wide In Vivo CNS Screening Identifies Genes that Modify CNS Neuronal Survival and mHTT Toxicity.

Unbiased in vivo genome-wide genetic screening is a powerful approach to elucidate new molecular mechanisms, but such screening has not been possible to perform in the mammalian central nervous system (CNS). Here, we report the results of the first genome-wide genetic screens in the CNS using both short hairpin RNA (shRNA) and CRISPR libraries. Our screens identify many classes of CNS neuronal essential genes and demonstrate that CNS neurons are particularly sensitive not only to perturbations to synaptic processes but also autophagy, proteostasis, mRNA processing, and mitochondrial function. These results reveal a molecular logic for the common implication of these pathways across multiple neurodegenerative diseases. To further identify disease-relevant genetic modifiers, we applied our screening approach to two mouse models of Huntington's disease (HD). Top mutant huntingtin toxicity modifier genes included several Nme genes and several genes involved in methylation-dependent chromatin silencing and dopamine signaling, results that reveal new HD therapeutic target pathways.

Genetic screens in isogenic mammalian cell lines without single cell cloning.

Isogenic pairs of cell lines, which differ by a single genetic modification, are powerful tools for understanding gene function. Generating such pairs of mammalian cells, however, is labor-intensive, time-consuming, and, in some cell types, essentially impossible. Here, we present an approach to create isogenic pairs of cells that avoids single cell cloning, and screen these pairs with genome-wide CRISPR-Cas9 libraries to generate genetic interaction maps. We query the anti-apoptotic genes BCL2L1 and MCL1, and the DNA damage repair gene PARP1, identifying both expected and uncharacterized buffering and synthetic lethal interactions. Additionally, we compare acute CRISPR-based knockout, single cell clones, and small-molecule inhibition. We observe that, while the approaches provide largely overlapping information, differences emerge, highlighting an important consideration when employing genetic screens to identify and characterize potential drug targets. We anticipate that this methodology will be broadly useful to comprehensively study gene function across many contexts.

Optimized libraries for CRISPR-Cas9 genetic screens with multiple modalities.

The creation of genome-wide libraries for CRISPR knockout (CRISPRko), interference (CRISPRi), and activation (CRISPRa) has enabled the systematic interrogation of gene function. Here, we show that our recently-described CRISPRko library (Brunello) is more effective than previously published libraries at distinguishing essential and non-essential genes, providing approximately the same perturbation-level performance improvement over GeCKO libraries as GeCKO provided over RNAi. Additionally, we present genome-wide libraries for CRISPRi (Dolcetto) and CRISPRa (Calabrese), and show in negative selection screens that Dolcetto, with fewer sgRNAs per gene, outperforms existing CRISPRi libraries and achieves comparable performance to CRISPRko in detecting essential genes. We also perform positive selection CRISPRa screens and demonstrate that Calabrese outperforms the SAM approach at identifying vemurafenib resistance genes. We further compare CRISPRa to genome-scale libraries of open reading frames (ORFs). Together, these libraries represent a suite of genome-wide tools to efficiently interrogate gene function with multiple modalities.

Uncoupling of sgRNAs from their associated barcodes during PCR amplification of combinatorial CRISPR screens.

Many implementations of pooled screens in mammalian cells rely on linking an element of interest to a barcode, with the latter subsequently quantitated by next generation sequencing. However, substantial uncoupling between these paired elements during lentiviral production has been reported, especially as the distance between elements increases. We detail that PCR amplification is another major source of uncoupling, and becomes more pronounced with increased amounts of DNA template molecules and PCR cycles. To lessen uncoupling in systems that use paired elements for detection, we recommend minimizing the distance between elements, using low and equal template DNA inputs for plasmid and genomic DNA during PCR, and minimizing the number of PCR cycles. We also present a vector design for conducting combinatorial CRISPR screens that enables accurate barcode-based detection with a single short sequencing read and minimal uncoupling.

Orthologous CRISPR-Cas9 enzymes for combinatorial genetic screens.

Combinatorial genetic screening using CRISPR-Cas9 is a useful approach to uncover redundant genes and to explore complex gene networks. However, current methods suffer from interference between the single-guide RNAs (sgRNAs) and from limited gene targeting activity. To increase the efficiency of combinatorial screening, we employ orthogonal Cas9 enzymes from Staphylococcus aureus and Streptococcus pyogenes. We used machine learning to establish S. aureus Cas9 sgRNA design rules and paired S. aureus Cas9 with S. pyogenes Cas9 to achieve dual targeting in a high fraction of cells. We also developed a lentiviral vector and cloning strategy to generate high-complexity pooled dual-knockout libraries to identify synthetic lethal and buffering gene pairs across multiple cell types, including MAPK pathway genes and apoptotic genes. Our orthologous approach also enabled a screen combining gene knockouts with transcriptional activation, which revealed genetic interactions with TP53. The "Big Papi" (paired aureus and pyogenes for interactions) approach described here will be widely applicable for the study of combinatorial phenotypes.

Creation of Novel Protein Variants with CRISPR/Cas9-Mediated Mutagenesis: Turning a Screening By-Product into a Discovery Tool.

CRISPR/Cas9 screening has proven to be a versatile tool for genomics research. Based on unexpected results from a genome-wide screen, we developed a CRISPR/Cas9-mediated approach to mutagenesis, exploiting the allelic diversity generated by error-prone non-homologous end-joining (NHEJ) to identify novel gain-of-function and drug resistant alleles of the MAPK signaling pathway genes MEK1 and BRAF. We define the parameters of a scalable technique to easily generate cell populations containing thousands of endogenous allelic variants to map gene functions. Further, these results highlight an unexpected but important phenomenon, that Cas9-induced gain-of-function alleles are an inherent by-product of normal Cas9 loss-of-function screens and should be investigated during analysis of data from large-scale positive selection screens.

Optimized sgRNA design to maximize activity and minimize off-target effects of CRISPR-Cas9.

CRISPR-Cas9-based genetic screens are a powerful new tool in biology. By simply altering the sequence of the single-guide RNA (sgRNA), one can reprogram Cas9 to target different sites in the genome with relative ease, but the on-target activity and off-target effects of individual sgRNAs can vary widely. Here, we use recently devised sgRNA design rules to create human and mouse genome-wide libraries, perform positive and negative selection screens and observe that the use of these rules produced improved results. Additionally, we profile the off-target activity of thousands of sgRNAs and develop a metric to predict off-target sites. We incorporate these findings from large-scale, empirical data to improve our computational design rules and create optimized sgRNA libraries that maximize on-target activity and minimize off-target effects to enable more effective and efficient genetic screens and genome engineering.

Rational design of highly active sgRNAs for CRISPR-Cas9-mediated gene inactivation.

Components of the prokaryotic clustered, regularly interspaced, short palindromic repeats (CRISPR) loci have recently been repurposed for use in mammalian cells. The CRISPR-associated (Cas)9 can be programmed with a single guide RNA (sgRNA) to generate site-specific DNA breaks, but there are few known rules governing on-target efficacy of this system. We created a pool of sgRNAs, tiling across all possible target sites of a panel of six endogenous mouse and three endogenous human genes and quantitatively assessed their ability to produce null alleles of their target gene by antibody staining and flow cytometry. We discovered sequence features that improved activity, including a further optimization of the protospacer-adjacent motif (PAM) of Streptococcus pyogenes Cas9. The results from 1,841 sgRNAs were used to construct a predictive model of sgRNA activity to improve sgRNA design for gene editing and genetic screens. We provide an online tool for the design of highly active sgRNAs for any gene of interest.